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ABSTRACT. Coastal areas, lagoons and river estu- 
aries along the Crimean peninsula, Sea of Azov and 
Russian Krasnodar region Black Sea coast have been 
searched for the presence of Harris mud crab [Rhithro- 
panopeus harrisii (Gould, 1841)] using snorkeling and 
SCUBA equipment. New populations were observed 
in 2 river estuaries (Vulan and Shap’sukho), Bugas 
liman adjacent to the Black Sea and it was widespread 
in the Taman Bay in the Sea of Azov. The crabs have 
been found in very distinct habitats from shallows of the 
Taman Bay with broad Zostera spp. meadows, to nar- 
row reed roots bands along the river estuaries. Similar 
habitats between the Black and Azov seas have shown 
complete absence of the crab in one and presence in the 
other respectively. Our results show that currently no 
standard method that exists, can be applied to all differ- 
ent types of habitats that we observed in this study. 
Nevertheless, the method of underwater transect census 
allowed us to perform quantitative analysis of crabs’ 
hiding behavior in the Taman Bay; smaller crabs of both 
sexes and most of ovigerous females prefer to hide in 
burrows. Size composition of crabs has been described 
for the Taman Bay and Black Sea estuary crab popula- 
tions. These studies lead us to suppose that there is a 
seasonal migration pattern of R. harrisii in the shallows 
of the Taman Bay and no such observations are appar- 
ent in the estuary habitat. While distribution of differ- 
ent age groups and migration patterns may be very 
different, the general life cycle of Harris mud crab in 
these populations appears to be similar. 


PE3IOME. IIpuo6pexusie yuacTKn, JlaryHbl M 39cTya- 
puu pek BAO Mobepexba Kpprma, A30BcKoro Mops H 


KpacHoylapcKoro moOepexba UepHoro Mops Opuin 06- 
CJI€OBAHBI C TOMOLI{bIO TOXBOAHOTO CHapsKeHHA C IC- 
JIO Haxon ena Rhithropanopeus harrisii (Gould, 1841). 
He w3BecTHbIe panee MOMyIAWMH OWI OOHapyKeHBI B 
actyapuax pek Bynas n Iancyxo n byra3ckom JiMMaHe 
Kpacuoylapcxoro kpas. KpaOpi oOvTaroT B CaMBIX pa3- 
JIM4HBIX MeCTOOOHTaHHAX OT MEJIKOBOJMÄĂ TamaHcKoro 
3asIMBa C NOJBOJHÞIMH JtyraMu Zostera spp. HO y3KHXx 
TOJIOC 3apocuel TPOCTHHKA B peuHbIX 9cTyapHAax. CxoL- 
HbIe MecTooOnTaHHa B UepHom H A30BCKOM MOpsX Xa- 
pakTepH30BasIHCb OTCYTCTBHeM BHMA B IIepBOM M IIpH- 
CyTCTBHeM BO BTOpoM. JlIMTebHbIe MOMyJIALHOHHEIC 
HCCJICNOBAHHA TpoBeyeHb! B TaMaHCKOM 3asIMBe C Ie- 
JIbIO pa3paOOTKH KOJIHYECTBEHHOTO MeTOa COopa Kpa- 
6os. Hanm pesyJIbTATSI MOKa3bIBaIOT, YTO TAKOTO EMH- 
HOrO MeTOJ{a, KOTOpÞIÑ MOXKHO ÕÞIJIO ObI TPEMJIOKHTE 
UIA BCEX HCCIeOBaHHbIX MECTOOÕHTAHHĂ, IIOKa He CY- 
INECTBYeT. TeM He mence, MCINOJIb3O0BAHHDIÑ METON TOJ- 
BOJHOTO TpaHCeKTHOTO yYeTa NO3BOJINJI pPOBECTH KO- 
JIMYCCTBCHHBIM AHaAJIM3 MCIIOIb30BaHHA KpaðaMH yKpBI- 
Tu: Kpa6sı Oosee MesKOorO pa3Mepa H OOJIbIIMHCTBO 
AMICHOCHBIX CaMOK TATOTEIOT K IIpe6bIBAHHIO B HOpKax 
H Apyrux yKphITHAX. [ia nonynannă Tamanckoro 3a- 
JIMBa H 9CTyYApHEB YeEPHOMOPCKOTO NobepexkbA ONCA- 
Ha pa3MepHaa CTpykTypa. ITO MCCJICMOBAHHE IO03BO0- 
JIET MpeIOJAraTb HAJIHYHE Ce30HHBIX MurpauHi R. 
harrisii Ha MmenkoBombax Tamanckoro 3amBa; HaOJIIO- 
JEHHA, NOLTBEpK LarolHe AHAJIOTHYHOE ABJICHHA Y Kpa- 
60B B 9CTyApHAX peK OTCYTCTBYIOT. B TO BpeMa Kak 
pacnpenenenne pasJIHYHbIX BO3PaCTHBIX Tpyill H HX 
MHTpaQHOHHaAA AKTHBHOCTBE B THX MOMMYIAWHAX OTJIH- 
yaeTca, OLHE YepThI KH3HCEHHOTO MKIA Kpaða Xap- 
puca NpegcCTaBJIAIOTCA CXOQHbIMH. 
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Introduction 


A small mud crab Rhithropanopeus harrisii (Gould, 
1841) or Harris mud crab originates from the eastern 
coast of American continent, from New Brunswick to 
northern Mexico [Williams, 1984]. In the native distribu- 
tion range it is generally associated with tidal estuaries 
where larvae of R. harrisii can maintain their distribution 
in estuarine areas by a complex mechanism of hatching 
synchronization with the tidal cycle and vertical migration 
superimposed to differences in velocity and directions of 
tidal flow [Forward, 2009 and references therein]. 

R. harrisii has invaded multiple basins around the 
world on almost all continents. In Europe (Fig. 1A), 
this mud crab species has been observed for the first 
time in 1874, in Zuiderzee, Netherlands [Maitland, 
1874; Buitendijk, Holthuis, 1949]. Since then it has 
spread along river estuaries and bays of the North 
Europe [Reznichenko, 1967; Van der Velde et al., 2000; 
Kerckhof, 2001; Olesen, Tendal, 2009], but only re- 
cently populations of this species were found in Lithua- 
nia [Baceviéius, Gasitinaité, 2008], Estonia [Kotta, 
Ojaveer, 2012], and Finland [Karhilahti, 2010; Fowler 
et al., 2013]. It has spread further south along the 
Atlantic coast of Europe [Eno et al., 1997; Mariscal et 
al., 1991; Cuesta et al., 1991; Projecto-Garcia et al., 
2010] and in to the Mediterranean [Mizzan, Zanella, 
1996; Noél, 2001; Ben Souissi et al., 2004]. 

In the Black Sea R. harrisii has been found in the 
Dnieper and Bug estuaries in 1937, much earlier than 
in the Mediterranean [Makarov, 1939]. It had then 
spread to several neighboring estuaries, and around 
1954-1957, Harris mud crab was reported from the 
open north-western Black Sea around Odessa [Rez- 
nichenko, 1967]. Simultaneously it was found through- 
out the Romanian and Bulgarian coasts [Zemiankow- 
ski, 1951; Marchand, Saudray, 1971]. Along the Crime- 
an coast of the Black Sea it has only been reported 
from the Sevastopol Bay [Shalovenkov, 2005]. In the 
Sea of Azov Harris mud crab was first observed in the 
Taganrog Bay in 1948 [Mordukhay-Boltovskoy, 1952]. 


By 1950s this crab became a common species through- 
out this sea [Reznichenko, 1967]. Consequently R. har- 
risii penetrated further eastward and became abundant 
in the Caspian Sea presumably through the Don—Volga 
Canal [Nebolsina, 1959; Reznichenko, 1967] and was 
then unintentionally introduced to the Aral Sea [An- 
dreev, Andreeva, 1988]. 

Surprisingly, so far, this species has not been report- 
ed in the published sources from the Taman Bay on the 
boundary of the Black and Azov seas, nor from the 
eastern Black Sea [except larvae found in some harbors, 
i.e. Selifonova, 2012]. The lack of data from such an 
extensive area in close vicinity to the old area of inva- 
sion limits precise mapping of Harris mud crab distribu- 
tion in the north-eastern Europe and restricts our under- 
standing of its ways of establishing new populations 
and the adaptive potential of this invasive species. 

Most studies that have looked at distribution and ecol- 
ogy of R. harrisii, have done so by collecting crabs using 
bottom dredges [i.e. Ryan, 1956; Turoboisky, 1973; He- 
gele-Drywa et al., 2014] and exposing crates at long time 
within natural substrate [i.e. Roche, Torchin, 2007; Fowl- 
er et al.,2013]. In the Black Sea these methods are diffi- 
cult to apply for rapid survey of an extensive area. Benthic 
substrates in the coastal zone are unsuited for dredging 
(i.e. rock, gravel or dense macrophytes beds) while long 
time setting of crates is impossible for most shores due to 
high exposure to wave impact. We therefore applied a 
methodology of in situ observations and underwater 
transect censuses, which in various versions is used in 
studies of commercial crab species [Miljutin, 2005]. 

The purpose of this study is the following: 1) to 
document occurrence of R. harrisii and its habitat pref- 
erence along the Crimean Peninsula, the boundary area 
between the Sea of Azov and the Black Sea and along 
the eastern coast of the Black Sea; 2) to test the appli- 
cation of direct transect census for studying population 
dynamics of this species; 3) to characterize local as- 
pects of sex and size composition of the sampled popu- 
lations; 4) to describe the behavior and habitat prefer- 
ences of the Harris mud crab in studied areas. 


Fig. 1. Distribution of Rhithropanopeus harrisii: A. Historical findings of Rhithropanopeus harrisii in Europe. For sources see text. 


Diamond — area of initial introduction in European waters in the XIX century; circles — areas where Rhithropanopeus harrisii crab has 
been found; Hashed rectangle — area of the present study. B. Black Sea and Sea of Azov regions where R. harrisii has been searched in the 
present study. Circles — spots searched, but no R. harrisii has been found; stars — R. harrisii has been found; dashed rectangles — data 
from these regions has been used for population size structure analysis in this study; grey region — R. harrisii has been known from this 
area previous to this study; square — recent records by other authors; 1 — Tarkhankut peninsula, 2 — Kerch, 3 — Cape Zheleznyi Rog, 4 — 
Anapa, 5 — Dyurso river, 6 — Yashamba river, 7 — Nechep’sukho river, 8 — Tuapse [Selifonova, 2012], 9 — Shepsi river, 10 — 
Matrosskaya river, 11 — Uch-Dere, 12 — Khosta. C. Study area in the Taman Bay, Sea of Azov. Numerated circles — stations of the 
benthic survey in 2008; empty circles — records of small juveniles of R. harrisii; filled circles — no crabs have been found; star — main 
area of yearly observations in transects in summers 2011-2014. 

Puc. 1. Pacupocrpanenue Rhithropanopeus harrisii B Espone: A. Uctopwueckne naxogkn BHA B Espone. MCTOYHHKH yka3aHbi B 
texcte. Pomé — paiion nmepBOHaAYAJIbHOTO BCeIeHHA B EBpPONCEÑCKHE BOMbI B XIX Beke: Kpyru — paŭňomsi, re R. harrisii Obi Halizeu 
BHOCICACTBUH; 3AITpHXOBAHHDIÄ IpPAMOYTONÞHHK — palioH WaHHoro nccnegoBanna. B. Mecra noncka Rhithropanopeus harrisii B 
UepHom H A30BCKOM MOpAX B XOJE MAHHOTO ucceqoBaHua. Kpyru — o6cIegoBannbie palioubl, rae kpa Xappuca He o6HapyXKeH; 
3Be3JO0YKH — paliOHbI HaxXoQOK Kpada Xappuca; 3aITpHXOBAHHAA OÕJACTb — aHHbIle H3 ITHX MECT OBLIM HCINMOJNÞ3OBAHDI JJA AHAJIH3A 
pa3MepHoli coctaBa kpaa Xappuica; palioubl c cepoli 3aIMBKOŘ — KpaO ŐbIJI H3BeCTeH H3 ƏTHX palioHOB panee; KBaypaT — y4acTOK 
HeJABHEÑ HAXOJKH Apyrumu aptopamu; | — Mbic Tapxankyr, 2 — Kepub, 3 — Mbic 3eneuniit Por, 4 — Anana, 5 — peka J[ropco, 6 — 
Pexa AniamOa, 7 — pexa Heuencyxo, 8 — Tyance [Selifonova, 2012], 9 — pexa Menca, 10 — pexa Martpoccxaa, 11 — Yu-Jlepe, 12 — 
Xocta. C. Paiion uccnegqoBanuii B TamMaHcKom 3amuBe A30BCKOTO Mops. IIponyMeposaHHbie kpyrH — OeuTocHbIe crann 2008 r. 
(cBeTIbIe KpyrH — kpaO ObuI OOHapyxXeH, TeMHBIe Kpyru — KpaO He OBL OOHapy2xKeH); 3Be3q04Ka — OCHOBHOe MeCTO TpaHCeKTHBIX 
yueTos KpaðoB B JeTHHM ce30H 2011-2014 rr. 
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Material and methods 


Study areas 


We have searched for the presence of Rhithropan- 
opeus harrisii in rivers and bays along the Black and 
Azov seas’ shores of the Crimean Peninsula and in 
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Krasnodar Region of Russia and the Taman Bay. (Fig. 
1B, Table 1). Only habitats where the Harris mud crab 
has been found are described further. For details of 
other spots see Table 1. 

Black Sea. The Bugas liman (Table 1; Fig. 1B) is a 
polyhaline (salinity 22%o) lagoon separated from the 
Black Sea by a narrow sandy strip and has a very 
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Table 1. List of studied localities in the Black Sea and the Sea of Azov. 
Taodmuua 1. Hepesenb oOcreqoBaHHbIx yY4ACTKOB B UepHom H A30BCKOM MOpsAX. 


R- Other crab 
Sea Region Spot, depth Position Year harrisii : Notes 
found species 
oaz z Pilumnus Closed bay; 

Black Sea | Crimea ee ol bay, AlS 5 E e i Jun 2012 No hirtellus mud, sand and 
shell 
Closed inlet; 

Sevastopol EAE mud, shell and 
A 3 °35'53.67", stones; 

Black Sea | Crimea oo river, E 33°36'26.98" Jun 2012 No Not recorded indications of 
significant 
eutrophication 

Black Sea | Crimea Sevastopol shore, Multiple places Jun 2012 No Eriphia Variety of 

0-10 m verrucosa, habitats 
ron , Liocarcinus Open shore: 

Black Sea | Crimea Tarkhankut X P O Jun 2012 No vernalis, P. pebble and 

peninsula, 0-5m | E 32°33'21.84 hirtellus rocky substrates 
There was 
flooding event 
in 2012, which 
dramatically 
rasaodas changed the 
Region Nechep’sukho N 44°14'34.77", banks of the 
Black Sea (Novomikhai- river 01.5 in E 38°50'12.82" Jul 2013 No Not recorded oe of 
lovey”) overgrowth and 
soft sediments 
have been 
flushed out; 
salinity 4%o 
Krasnodar Shap’sukho N 44°18'13.40", | Aug 2012 
Black Sea | Region river** , 0-2 m E 38°44'38.40" Jul 2013 , Yes Not recorded 
(Lermontovo *) 
Krasnodar 
Black Sea ae Mee i as : pa e ie E Yes Not recorded Ha text for 
: etails 
Osipovka *) 
Krasnodar 
Black Sea ee Hee liman, ae S ‘ Bes á Jun 2014 Yes Not recorded 
skaya*) 
Open shore 
areas dominated 
by sandy 
Brachynotus beaches and 
sexdentatus** | extensive 
*, Carcinus shallow sandy 
Krasnodar Aug 2012 aestuarii, L. areas. At 
Black Sea Region Bugas peninsula N 45° 3'53.52", Jul 2013 : No vernalis, P. approximately 
(Blagoveshchen- | sea shore, 0-5m | E 37° 2'21.07" J 201. 4 hirtellus, 2-5 meter depth 
une A : 
skaya*) Macropodia there is a sandy 
czernjawskii, bottom with 
Xantho crushed shell 
poressa, debris and 
patches of 
filamentous 
algae 
Open sandy and 
Krasnodar 625% ss ebble beach, 
Black Sea | Region 0-4 m N ar 3335:07" , Aug 2012 No Not recorded indications of 
E 37°18'48.14 
(Anapa *) strong 
eutrophication 
Clay cliffs and 
E. verrucosa, benches with 
Krasnodar Near cape N 45° 6'50.56" L. vernalis, reefs (with 
Black Sea | Region Zheleznii Rog, E 36°45'0 93" á Aug 2012 No M. bores of a rock- 
(Volna *) 0-3 m i czernjawskii, boring mollusk 
P. hirtellus Pholas 
dactylus) 
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Table 1 (contituing) 
Tadsmua | (mpozomxKeHne) 


Krasnodar earn ‘i 
Black Sea | Region 0-2 m Aiea: k Aug 2012 No 
(Khosta *) f 
Krasnodar 020! " 
Black Sea | Region 0-2 m anes es ae Aug 2012 No 
(Uch-Dere *) ` Open rocky 
Krasnodar Matrosskaya N 43°48'0.05" Hore 
Black Sea | Region : y: SBAA AA Aug 2012 No intermitted by 
: river**, 0-2 m E 39°27'12.10 E. verrucosa, 
(Golovinka *) pebble beaches; 
Pachygrapsus $ : 
Krasnodar ki P ‘i rivers with 
; Shepsi river**, N 44° 1'59.71", marmoratus 
Black Sea | Region 02 E 39° 8'30.85" Aug 2012 No poorly 
(Shepsi *) ot i developed 
Krasnodar ; err P estuaries 
Black Sea | Region ya (e SAKE en ý Aug 2012 No 
(Gelendzjik *) ij ` 
Krasnodar s 
: Dyurso river**, N 44°40'42.00", 
Black Sea | Region 02m E 37°33'34.76" Aug 2012 No 
(Dyurso *) 
Bad weather 
Sea of . $ N 45°19'16.15", i 
‘hee Crimea Kerch E 36°28'52.40" Jun 2012 No Not recorded conditions 
during survey 
Sea of ‘i # N 45°26'37.17", See text for 
Ao Crimea Zolotoe*, 0-2 m E 36° 4'31.66" Jun 2012 Yes Not recorded details 
3 Shallow area 
Sea of Krasnodar Beysug Liman, N 46° 9'13.09", : : 
Azov Region oim E 38°28'1.50" Aug 2012 No Not recorded | with extensive 
reed coverage 
March 
2008 
: Jul 2008 
Multiple spots on 
Taman Krasnodar 0-4m ihe-codst and Jul 2011 Yes P. hirtellus 
Bay Region insidethe ba Aug 2012 
ue y Jul 2013 
Apr 2014 
Jun 2014 
Jul 2008 
Jul 2011 
Taman Krasnodar Coastal shallows | N 45°16'43,09", | Aug 2012 , 
Region in inner part of dimi ‘i Yes P. hirtellus 
Bay (Sennoi *) bay, 0-2 m E 36°57'45.28 Jul 2013 
> Apr 2014 See text for 
Jun 2014 details 
Taman Krasnodar Rubanov N 45°19'28.77", ; 
Bay Region Peninsula, 0-2 m | E 36°47'9.67" nents xes G S 
Krasnodar arai ‘i 
Taman Region 0-2 m N ay kai oe j Jul 2013 Yes P. hirtellus 
Bay E 36°44'31.53 
(Taman *) 
Taman Krasnodar : N 45°14'40.91", Jul 2009, 
Bay Region Tuzla spit, 0-3 m E 36°35'28.98" Aug 2012 Yes Not recorded 
Krasnodar 
Taman Region Inner part of bay, | N 45°19'43.31", 
Bay (Volna Revo- 0-2 m E 36°55'50.08" | 112013 Yes [-Norrecorged 
lutsii *) 


Legend: * Villages/Towns, ** Both river estuary and nearby sea shore has been searched. *** first record of the species in this region 
of the Black Sea. See Figures 1 B,C 


narrow canal that connects these two water bodies. It is 
very shallow and has soft bottom sediments covered by 
a thick bacterial mat. There are sporadic patches of 
macrophytes and large numbers of /dotea balthica and 
Cerastoderma glaucum. 

Estuaries of major rivers along the shore of the 
Black Sea in Krasnodar region, Russia have been in- 
spected. Most of these rivers are subject to large flood- 
ing events during the rainy season, and very low water 


levels in summer and near 0%o salinity. Rivers 
Shap’sukho near Lermontovo village and Vulan near 
Arkhipo-Osipovka village (Table 1; Fig. 1B) are sub- 
ject to less flooding during rain. Estuaries of these 
rivers unlike most others have muddy shores and a 
lower flow speed which allows salty water to enter 
from the sea raising their salinity levels up to 9%o. 
These rivers have reed growth along its shore and 
islands. 
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Sea of Azov. An open shore near village Zolotoe 
on the Sea of Azov side of the Crimean Peninsula 
(Table 1; Fig. 1B) is subject to high wave impact. 
Sediments were mostly crushed shell debris and stones, 
with sparse macrophyte vegetation (Fig. 2A). 

Taman Bay. The Taman Bay is formally included 
in the Sea of Azov, but it has a distinct hydrological 
regime. It became a regular area of study (between 
2008 and 2014) after the black oil spill in the Kerch 
Strait in November 2007 [see Kolyuchkina et al., 2008, 
2009, 2012] (Fig. 1C). This bay is a large (about 350 
km’), shallow (maximum depth does not exceed 5 m) 
embayment which opens to the Kerch Strait in the 
south-western part of the Sea of Azov. In the western 
part it is separated from the Kerch Strait by Chushka 
spit and in the south by Tuzla spit which has been 
artificially extended and strengthened in mid-2000s. 
The bay lacks freshwater continental discharge; its sea 
level and hydrological regime are strongly influenced 
by dominating winds. Temperature may rise up to 27— 
28°C (and even up to 32°C as we have observed in 
shallow areas in 2009) in summer and fall, to a freezing 
point in winter. Salinity in the bay usually ranges be- 
tween 14.8 and 15.5%o, which is higher than the salini- 
ty of the southern Sea of Azov (10-11%) [Zalogin, 
Kosarev, 1999; Golovkina, Nabozhenko, 2012]. Our 
measurements indicate that it may reach up to 18%o, if 
winds drive the Black Sea water inside the bay. A 
peculiar characteristic of the Taman Bay is the impor- 
tance of macrophytes in the ecosystem and large shal- 
low water areas. 

Shallow zone near Sennoi village is very broad. 
Depth of 2 m is approximately 200 meters away from 
the shore (Fig. 1C). We have chosen this area for our 
study due to an easy access and its shallowness. Ap- 
proximately 100 m away from the shore it is still an 
arm’s length depth. This allowed us to collect crabs 
without any special equipment, just with a mask, to 
search thoroughly through the sediments as well as to 
observe crabs without scaring them with our move- 
ments. During weak wind conditions the water is clear 
and the majority of crabs can be visually spotted. In 
colder seasons strong winds can blow away large 
amounts of water resulting in open sediments of up to 
50 m from the shore [pers. comm. local community, 
pers. observations by AZ]. Sediments in this zone are 
silty sands, mostly covered by macrophytes (Potamo- 
geton sp., Zostera noltei, Cystoseira sp., green fila- 
mentous and brown algae) with patches of open sandy 
sediments. The top oxygenated layer is about 3 cm 
thick, below a dark, anoxic layer is present. Infaunal 
communities are dominated by cockle (Cerastoderma 
glaucum), another bivalve Loripes lucinalis and a gas- 
tropod Bittium reticulatum. For the epifauna, patches 
of blue mussels (Mytilus galloprovincialis), barnacles 
(Balanus improvisus) and sea anemones (Actiniaria, 
yet unidentified and probably an alien species; N. San- 
amyan, pers. comm.) on available substrates are char- 
acteristic. Active benthic macrofauna mainly consists 
of Gammarus aequicauda, Idotea balthica, R. harrisii, 


Palaemon adspersus; among benthic fishes gobies (Go- 
bius spp.) and Sygnathus nigrolineatus are common. 
Deeper than approximately 2 m, Z. noltei is substituted 
by eelagras meadows of Zostera marina. Deeper than 
3.54 m seabed is dominated by muddy sediments with 
thick bacterial mats, rare patches of macrophyte vege- 
tation, sea anemones and other sessile macrofauna sit- 
ting on dead cockle shells (underwater observations by 
the authors). 


Observations and collecting methods 


For short site visits along the Black and Azov seas’ 
shores we used a standardized observation scheme. The 
non-estuary shore sites were usually surveyed within a 
day or two using SCUBA diving or snorkeling equip- 
ment from the watermark to the depth of 2—5 m (depend- 
ing on the bottom topography). All underwater excur- 
sions were done during daytime and only if there was 
appropriate visibility. During the excursion we inspect- 
ed all characteristic biotopes, carefully examining those 
where crabs could be hiding, i.e. stones, holes, crevices, 
macrophyte beds, fouling communities, artificial sub- 
strates and shore protection constructions. 

In estuaries, where low visibility precluded under- 
water observation, we inspected shores, stones, litter, 
artificial substrates and reed plants using hand net 
and spade. 

In the Taman Bay we used data of macrobenthic 
survey (conducted between 7 and 29 July 2008) to 
reveal distribution of juvenile Harris mud crabs. Sam- 
ples have been taken using SCUBA diving equipment 
with manual tube sampler (0.0095 m°); 3 to 5 casts at 
each station. (Fig. 1C) Samples have been washed 
through a 0.5 mm mesh, fixed with buffered 4% form- 
aldehyde and later sorted in the laboratory. 

Transect censuses in the Taman Bay have been 
conducted near village Sennoi. R. harrisii has been 
sampled in July 2011, August 2012, July 2013 and 
June 2014. Crabs have been collected from 100 m? 
rectangular transect perpendicular to the shore. The 
first author walked 100 average steps away from the 
shore and swam back to the shore observing and col- 
lecting all crabs that were within 1 m reach around her. 
The maximum depth was at an arm’s length (approxi- 
mately 1 m). The last half a meter, near the shore, was 
often omitted due to low visibility. The researcher, 
using snorkeling equipment and a wet suit, spent up to 
an hour passing one transect, depending on the quanti- 
ty of crabs observed. All crabs that were not in burrows 
were placed in a separate sampling bag to those that 
were found in burrows. All large double or multiple 
burrows were poked with finger and crabs emerged 
from the opposite entrance. Roots of macrophytes were 
poked through with fingers in search for pockets where 
females might hide. In areas with substantial Zostera 
noltei coverage the grass has been “brushed” by hands. 
All empty bivalve shells and stones were overturned. 
Specimens from most sampling efforts in 2011 have 
been measured with a ruler on the shore and released. 
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Pe i j 
Fig. 2. Examples of habitats where Rhithropanopeus harrisii has been observed: A — sandy/shell debris habitat near the village 
Zolotoe, Sea of Azov. Crabs have been found hiding under some stones that were scattered on the bottom, with high density of gobies fish 
presence. B, C — burrows in the Zostera noltei medow near village Sennoi, Taman Bay, 2011. D — muddy sediments with thick bacterial 
mat in Bugas liman, near village Blagoveshchenskaya. E — patchy cover of macriphytes in the shallows near village Sennoi, Taman Bay in 
2013. 

Puc. 2. Ipnmeps MectooOutanuit Rhithropanopeus harrisii: A — necok n pakyma y moc. 30m0Toe, AsoBckoe Mope. Kpa6sı bum 
oOOHapy2KeHE! MOJ OTACIbHBIMH KAMHAMH; IIPHCYTCTBOBAsIO ÕOJBINOE KOJIHYECTBO ÕbIYKOB. B, C — HopKH B sapocax Zostera noltei y noc. 
Cennoŭ, TamaHckuit 3auB, B 2011 r. D — umuctoe HO c 6akTepHabHOÑ MJICHKOÑ B byrackoM JIHMAaHeE, y cTaHuubl baaroBenlenckas. 
E — M03aHYHbIÑ MOKpOB MakpoþHTOoB y moc. Cennoŭ, Tamanckuii 3am, 2013 r. 


In other years most of crabs have been fixed in 10% has been made due to high quantities of filamentous 
formaldehyde and transported to the laboratory. Crabs algae that restricted visibility. 

collected in August 2012 were not used in this study In river estuaries crabs have been collected by cut- 
since no successful completion of a standard transect ting off parts of “root walls” created by reeds along the 
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shore. Chunks of these roots have been detached with a 
spade directly into a bag, washed through and all crabs 
were picked out by hand. 


Analysis 


All crabs have been sexed, carapace width (CW) 
has been measured with accuracy 0.2 mm using a cali- 
per or a ruler (in 2011). For the analysis we used CW 
sizes of crabs collected by a standard method from the 
Taman Bay, near Sennoi village, juveniles found in the 
benthic 2008 Taman Bay samples and from river Vu- 
lan in Arkhipo-Osipovka, where the sampling method 
is likely to represent local population structure with 
low sampling error. To analyze population size struc- 
ture, a mixing model analysis of CW sizes has been 
performed using PAST software package [version 3.0; 
Hammer, 2013]. The Akaike information criterion 
[Akaike, 1974] and the log likelihood parameter have 
been used to identify best fit model. The mixture model 
analysis has been re-run 50 times. All ANOVA tests 
have been performed using R software package [R 
Core Team, 2013]. 


Results 


Rhithropanopeus harrisii occurrence 


Black Sea 

Crimean Peninsula. Site visits along the western 
coast of Crimea did not result in observations of Rhithro- 
panopeus harrisii (Table 1). 

The Black Sea coast in the Krasnodar region. 
No indications of R. harrisii presence in the open shore 
areas of Khosta, Uch-Dere, Golovinka, Shepsi, Ge- 
lendzhik, Dyurso, Anapa, Blagoveschenskaya, and Vol- 
na have been observed (Table 1; Fig. 1B). 

In the Bugas liman (Table 1; Fig. 1B), R. harrisii 
were sitting on macrophyte beds, hid in the folds of 
bacterial mat and rarely in burrows which were scarce 
and mostly populated by gobies (Fig. 2D). As men- 
tioned above, no R. harrisii have been observed in 
close vicinity to the liman on open shore area, although 
salinity in the liman was higher (22.5%) than that of 
the open sea (17%o), or the Taman Bay (14.8—15.5%o). 

Rivers 

Harris mud crabs have been found on an island in 
River Shap’sukho estuary (Table 1; Fig. 1B) over- 
grown by reed. Due to low visibility, we were only able 
to observe the very top layer of the island and no 
burrows have been observed. Crabs have been seen 
feeding on the overgrowth of submerged reeds. They 
were dark colored and hid at first sight of any move- 
ment. Only very small amount of specimens have been 
collected in August 2012 and no specimens have been 
caught in July 2013 although they have been spotted. 

In the estuary of River Vulan (Table 1; Fig. 1B) 
Harris mud crabs have been found in rare shallow areas 
under stones. More commonly they have been collect- 
ed from the roots of reeds along the shore. These roots 


created an abrupt wall with deep network of crevices 
where the crabs resided. All crabs were dark colored 
and hid at any sign of the collector’s movement. Not 
much visual observations could be made due to very 
low visibility. Bottom sediments along these walls were 
very soft and bubbles of gas came out when pressed. 

Taman Bay. During the underwater sampling and 
excursions, Harris mud crab has been spotted in most 
parts of the Taman Bay, i.e. in the coastal area of Tuzla 
spit, near Taman town, in the central deepest part, in 
the inner part near Sennoi, Volna Revolutsii and 
Garkusha villages at depth between 0 and 4 m (Table 
1; Fig. 1C). It occurred in practically all kinds of habi- 
tats known there: sandy bottom with different amount 
of mud and shells and mosaic macrophyte vegetation 
and philamentous algae, in seagrass meadows, on clay 
reefs and patches with a labyrinth of crevices, in mus- 
sel beds, on muddy bottom with cockle shell and sparse 
macrophytes, and on stone blocks of the antique jetty 
exposed during underwater archaeological excavations 
by the Fanagorian Archaeological Expedition in Sen- 
noi (Fig. 2B, C, E). Of particular interest are the early 
spring (8.03.2008) observations made in front of Sen- 
noi where relatively numerous crabs have been ob- 
served at depth 3.2 m among cockles and sparse eel- 
grass on the grey-brown mud with a thin (about 1 cm) 
upper light oxygenated layer and a thick dark anoxic 
layer under-laid by dense shell deposits. The only area 
where we have not spotted R. harrisii is the very shal- 
low (usually less than 1.5 m) north-western part of the 
Bay, in so called Dinskoi Bay and along the Chushka 
spit where extensive areas are occupied by reed [Spiri- 
donov et al., 2016]. 

Small juvenile Harris mud crabs were practically 
absent in our hand-picked samples, while the benthic 
survey conducted in the Taman Bay in July 2008 in 
deeper areas (Fig. 1C) revealed them to be present at 
30% of stations. They were however absent in the very 
shallow north-western part of the bay. In the eastern 
part of the Taman Bay their frequency of occurrence 
was 43%. Small sizes of most crabs (see below) and 
the presence of megalopa larvae and early crab stages 
indicate that they recently passed metamorphosis and 
settled on the substrata. In most spots, where juvenile 
crab were present, sediments were muddy (alevropelite 
fraction exceeded 50%) except stations 23 and 26 (Fig. 
1C) where dense shell substrate has been observed. In 
half of these stations juvenile crabs settled in dense 
eelgrass meadows, in other cases seagrass was patchy 
but seabed was covered with green or brown algae. 
Macrobenthic biomass in all stations where crabs were 
present was high (217-1573 g per m°). Particularly 
high values have been observed for the communities 
dominated by Mytilaster lineatus, C. glaucum and An- 
adara kagoshimensis [Kolyuchkina et al., 2013]. At 
station 26, adult crabs along with numerous sea anemo- 
nes have also been observed. Rough calculation indi- 
cate that the number of just settled mud crabs in eel- 
grass belt may be as high as 125-325 individuals per 
m’, but these figures have to be taken with caution 
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Table 2. General information on standard transects for censuses of Rhitropanopeus harrisii in the Taman Bay, Sea of Azov 


in different years. 


TaOmmua 2. Omas nuþopmarnma mo CTaHJapTHBIM TpaHcekTaM, Ha KOTOPBIX MpoBogzMIMYb yeri Rhitropanopeus harrisii 


B pa3Hble roga B TaMaHCKoM 3aJIHBe A30BCKOrO MOp4. 


GPS Year Salt Temp Not Nor Crab per Notes 
Cc transects crabs m 
10-18 Jul 2011 16,5 8 1096 1,440 
no visibility, high 
Aug 2012* 0 20 content of 
filamentous algae 
N 5°16'43.09" = oe 2 and 
E 36°57'45.28" | 4-12 Jul 2013** 19,5 24 7 64 0,091 i a 
3 (no crabs found) 
meter depths 
20-25 April 
2014* 17 19,5 6 20 0,033 
19-23 Jun 2014 17 23,5 4 21 0,053 


Legend: * Not used in this paper, ** 5 standard and 2 parallel 


depth parallel transect). 


owing to sampling constraints. In July 2011 we occa- 
sionally observed very small R. harrisii in eelgrass 
meadows at depth of about 2.5 m in front of Sennoi. A 
benthic survey of 15 stations, using the same methods 
and having even broader coverage than the survey of 
2008, did not reveal presence of recently settled crabs 
in the Taman Bay in July 2013. 

Pilumnus hirtellus (Linnaeus, 1761) is the only oth- 
er brachyuran species found in the Taman Bay (Table 
1). Near Taman town and in the vicinity of Rubanov 
Peninsula where the sediments consist of open clay 
patches with a labyrinth of crevices R. harrisii shares 
these crevices with P. hirtellus, where they coexist side 
by side. 

The collections of the Zoological Institute of the 
Russian Academy of Sciences in St. Petersburg indi- 
cate an occurrence of Brachynotus sexdentatus (Risso, 
1827) in the Bay prior to the introduction of R. harrisii 
(catalogue number 38281). This species was not found 
during our 5-years studies in the Taman Bay in differ- 
ent seasons. 


Transect observations and censuses in shallow 
coastal flats of the Taman Bay 


Crabs’ habitat and habits 

Specimens of Rhithropanopeus harrisii in the shal- 
low water near village Sennoi are usually lightly col- 
ored (sandy yellow) and are often found in burrows as 
well as in the open. We have never observed this crab 
in the process of digging a burrow, thus we cannot 
assume they are made by the crab. Male or some large 
females of R. harrisii are often found in large burrows 
with 2 or multiple entrances (Figs 2B, C and 3). This is 
very similar to the burrows occupied by gobies. Rarely, 
they have been observed sharing a burrow with multi- 
ple entrances with the fish, but never with each other. 

Many females lead a cryptic life and in 2011 they 
were often found in small pockets underneath macro- 


transects are used in this paper’s analysis (no crabs found at 3 meter 


phyte roots, with no apparent entrances. These pockets 
were only located by pocking through the roots with 
fingers. None of these pockets were found in subse- 
quent years or in any other location. When staying in 
the open, most of large males were assuming a defense 
position with open claws raised towards the observer. 
Smaller crabs frequently flee to their nearby burrows 
or burry themselves in rare sandy patches, which could 
be spotted easily due to darker colored sand. Often 
crabs took no notice of the researcher and continued 
their feeding. During the observations the sampler hardly 
moved and took a long time to search the area. 

Crab counts 

After several trials in 2011 we have successfully 
completed 8 transect censuses: 3 at 6:30 in the morning 
(1.5 crabs per m°), 2 at noon (1.37 crabs per m°) and 3 
at 6:30 in the evening (1.4 crabs per m°), on average 
1.44 Harris mud crabs per m’. 

Overall we have collected 1096 living crabs: 434 
(40%) in the morning, 264 (24%) at noon and 398 
(36%) in the evening. More than half (about 60% of the 
total count) crabs have been found in burrows. Sex 
ratio was biased towards males 78%, 13% non-oviger- 
ous females and 10% were ovigerous females. The 
breakdown of numbers and percentages of crabs found 
in July 2011 are summarized in Tables 2 and 3. 


Table 3. Number of crabs found at different time of day on 
different standard transects in July 2011, in the Taman Bay. 
TaOmuua 3. Konmuectso Kpa6os oOHapy2KeHHBIX B pa3Hoe 
BpeMA CYTOK Ha CTaHapTHbIX TpaHceKTax B MONE 2011 r. 
B TaMaHCKOM 3aJIMBe. 


ome . the - = - Tol 
Morning 136 187 111 434 
Noon 121 143 264 
Evening 102 175 121 398 
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Table 4. Number and percentage of male, non-ovigerous (non-ov.) and ovigerous (ov.) females, found at different time of 
day on standard transects in July 2011, in the Taman Bay (sum in a row equals to 100%). 

Tamua 4. KomwuectsBo u gons (%) caMijoB, caMok 6e3 (non-ov.) H c HKpoli (ov.) Ha meonogax, OOHAPYKCHHBIX B 
pa3HOoe BpeMA CYTOK Ha CTaHapTHbIX TpaHceKTax B HONE 2011 r., B TamaHcKom 3as1mBe (CyMMa B cTpoke 100%). 


N of Males Non-ov. females Ov. females 

samples n % n % n % 
Morning 3 327 75% 61 14% 46 11% 
Noon 2 215 81% 30 11% 19 7% 
Evening 3 308 77% 47 12% 43 11% 

Total n 8 850 -— 138 -— 108 -— 


One-tailed ANOVA test has been performed to com- 
pare sizes of crabs found at particular transect census- 
es. There was a significant difference (df=7, p <2e-16), 
but a post Hoc Turkey test has revealed that there is 
only one transect which is an outlier. Crabs from this 
transect were measured in the lab using calipers where- 
as a ruler has been used to measure the CW of crabs 
from the 7 transects in the field. Due to this, the mea- 
surement error is not comparable. Without this sample 
all transect samples revealed no significant differences 
in CW of collected crabs (df=6, p=0.156); this indicat- 
ed a comparability of sampling effort between particu- 
lar censuses. Thus, for further analysis using size of 
specimens, the sample measure by the caliper was ex- 
cluded from the data (resulting number of observations 
909). 

In July 2013 no females have been found in the 
vicinity of Sennoi and Volna Revolyutsii villages and 
the Rubanov peninsula. Near the shore, close to Taman 
town, no females have been observed either, but fur- 
ther away, at approximately 4.5 m depth, several fe- 
males have been collected. Due to this finding 3 more 
100 m transects (at 1.5, 2 and 3 m depths) using SCU- 
BA instead of snorkeling have been made. These 
transects were parallel to the shore where the standard 
perpendicular transects have been sampled near Sen- 
noi village. Crabs were found only at 1.5 and 2 m. 


Once again, no females have been found at these depths. 
Overall 5 standard transects yielded 34 crabs, on aver- 
age 0,068 crabs per m?. This is a much lower figure 
when compared to the 2011 data. 

Similarly to July 2013, in June 2014 all crabs found 
near Sennoi were males except for one large female 
(CW 15.3 mm). During sampling period the weather 
was stormy with high waves, which destabilized bot- 
tom sediments in shallow waters. Only 4 samples have 
been successfully obtained which yielded 21 crabs 
(0.0525 crabs per m°). Thus summer abundance of 
mud crabs in coastal shallows near Sennoi decreased 
by two orders of magnitudes between the years 2011 
and 2013-2014 (Table 4). 

Hiding pattern of crabs with regard to sex, size and 
time of day 

The most abundant samples of 2011 make it possi- 
ble to analyze the crabs’ pattern of hiding in burrows. 
The fraction of crabs found in and out of burrows 
varied among transects, the time of the day and be- 
tween the males and females but the proportion of 
crabs hiding in burrows and other types of shelter i.e. 
pockets underneath macrophyte roots (see above) was 
always significant, 41 to 78%, 58.4+6.08% in average 
(Tables 5, 6). 

To reveal the differences in size of crabs found in 
and out of burrows and other shelters at different time 


Table 5. Number and percentage of male, non-ovigerous (Non-ov.) and ovigerous (Ov.) female of Rhithropanopeus har- 
risii found in (In) and out (Out) of burrows at different time of day on standard transects in July 2011, in the Taman Bay 


(sum in a column equals to 100%). 


Tamuna 5. Kommuectso H oma (%) caMNOB, caMoK 6e3 (Non-ov.) n c HKpoli Ha mueonoyax (Ov.) B mpo6ax Rhithropan- 
opeus harrisii, coOpaHHsbix B HopKax (In) n BHe HOpoK (Out) B pa3Hoe BpeMA CYTOK Ha CTaHapTHbIX TpaHCeKTAaXx B HIOJIC 


2011 r., B TamaHckom 3amBe (cyMMa no croxõuyam 100%). 


Males Non-ov. females Ov. females 

N 43 

n 35 
In 

% 81% 

n 8 
Out 

% 19% 


Mor. — Morning, Even. — Evening. N — total number of crabs in particular group. 
Mor . — ytpo, Even. — Beyep, Ov. — camxu C uKpoii. N — o6mee kommuecTBO Kpa0os B OlpeyeeHHO rpynne. 
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A 


Fig. 3. Schemes of burrows where Rhithropanopeus harrisii can be found in the Taman Bay, Sea of Azov. Drawing by A.K. Zalota 
based on personal observations: A — burrow with multiple entrances, B — burrow with two entrances, C — burrow under roots; in this 


type of shelter most of females have been found in 2011. 


Puc. 3. CxeMbI THIIOB HOpPOK H JPyrux yőeKHN, B KOTOPBIX Öb OOHapyxeH Rhithropanopeus harrisii B Tamanckom 3aMBe. PucyHoK 
A.K. 3anota Ha OcHOBe JIM4HBIX HaOmoqeHHit: A — HOpKa C HeCKObKHMM BXOMAMH, B — HopkKa C JByMa Bxogamu, C — HopkKa B 
KOpHeBHNJAX MaKpOUTOB; B yŐexumax TAKOTO THNAA OOHApy2KeHO OOJIIMHCTBO CAMOK B 2011 r. 


of day we have performed a two way ANOVA test 
(factor A — in vs out of burrows, factor B — samples). 
The results are summarized in Table 7 and Figs 4 and 5. 

The sizes of crabs do not differ significantly in 
most cases, except for the non-ovigerous females in the 
morning. From the box plot (Fig. 5) it is possible to 
conclude that in the morning smaller females are hid- 
ing in burrows whereas larger are out and are possibly 
foraging. There are no statistically meaningful results 
for ovigerous females since no specimens have been 
found out of burrows and only 3 specimens were found 
hiding in the morning. Similar circumstances were ob- 
served at noon, thus the ANOVA results for ovigerous 
females at noon should be treated with caution. 

One way ANOVA test has been used to analyze the 
differences in sizes of crabs collected between morning, 
noon and evening for each sex in and out of burrows. 
Similarly due to low sample size of ovigerous females 
found out of burrows no analysis has been done. Only 
ovigerous females in burrows have revealed significant 
difference of CW at different time of day (p=0.011), 
although based on the box plot (Fig. 5) we would expect 
to find such differences in all females. Post hoc analysis 
has revealed significant difference between the morning 
and the evening samples (p=0.0076). 

Mixing model analysis has indicated that the male 
population in burrows comprises of 5 size groups, while 
only two most numerous groups (corresponding to 
groups 3 and 4) are present in the males collected out 
of burrows (Table 8; Figs 6A, B). The smallest groups 
with mean CW 7.6 mm and 9.6 mm and the largest 


group (mean CW 18.6 mm) have not been recognized 
in males out of burrows (Table 8) which were however 
1.56 times less numerous (Table 7). Although the num- 
ber of non-ovigerous females in burrows was less than 
the number out of burrows, three groups have been 
recognized in this sample and only most numerous 
group 2 was present in the non-ovigerous females col- 
lected out of shelters (Table 8; Figs 6C, D). A very 
similar situation with the exception of the largest group 
absence was found in the ovigerous females (Table 8; 
Figs 6E, F). 

Overall from Tables 5 to 8 and Figs. 4—6 it is clear 
that the majority of males hide during the day (except 
in the morning similar numbers have been observed in 
and out of burrows). There are no major differences in 
sizes except the smallest males were more common in 
burrows. Non-ovigerous females preferred to be out 
during the day and became hidden towards the evening. 
There is a significant difference in size of females in 
the morning, where larger females were out of shelters 
and smaller ones were hiding in the burrows. Oviger- 
ous females generally had cryptic habit during day 
light and few of them were found outside during our 
study. 


Aspects of size composition 


Quantitative sampling in shallow coastal flats in 
Taman Bay 

For the largest sample collected near Sennoi in 2011 
the Mixture Model Analysis has shown three size groups 
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Fig. 4. Box plot graphs of sizes of crabs found in or out of burrows grouped by sex and time of day, collected on standard transects in 


July 2011, in the Taman Bay. 


Puc. 4. J[uarpamMbI pa3Maxa pa3MepoB Kpa0oB, HaliJjeHHbIX B H BHe HOpOK, (pa3eIeHHBIX MO NOY H BpeMeHH CyTOK; cOop Ha 


CTaHJapTHBIX TpaHceKTax B HONE 2011 r. B TamaHcKoM 3asIMBe. 


of males with mean CW 12.4, 15.6 and 18.5 respective- 
ly, with the largest group constituting a small fraction of 
the sample (Table 9; Fig. 7A). Two additional small size 
groups were also identified in the subsample collected 
from burrows (Table 8; Fig. 6A, B). Females on the 
other hand had only two size groups (mean CW 9.2 and 
12.1, where the second group in females corresponded 
to the group with mean CW 12.4 mm in males and the 
first group corresponded to the additional small group 
that was absent in the total male sample but present in 
the subsample collected from burrows). 

Both ovigerous and non ovigerous female popula- 
tions have been shown to comprise of the same two 


groups with characteristics similar to the overall fe- 
male sample, the smaller group was proportionally less 
represented (Table 9; Figs 7B—D). 

In July 2013 the sample was much smaller com- 
pared to 2011 and contained only males. For Mixture 
Model Analysis we used crabs that have been collected 
at the 5 standard and 2 depth transects. This analysis 
has produced 2 size groups with mean CW of 11.8 and 
18.7 mm with a clear domination of a larger group 
(Fig. 8B). These two groups corresponded respectively 
to the smallest and the largest group of males in 2011 
while anything similar to the intermediate group of 
2011 was not detected. 
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Males In crabs' CW vs time 


Female In crabs’ CW vs time 


Fov In crabs’ CW vs time 
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Fig. 5. Box plot graphs of CW sizes of crabs collected at different time of day grouped by sex and in or out of burrows, on standard 


transects in July 2011, in the Taman Bay. 


Puc. 5. J[marpammsi pa3Maxa pa3MepoB Kpa0os, HaliJIeHHBbIX B pa3HOe BpEMA CYTOK (pa3eJICHHBIX IMO NONY H HAXOKICHHIO B HJIH BHE 
Hop); COop Ha cTaHJapTHBIX TpaHceKTax B HIONE 2011 r. B TamaHcKoM 3aJnBe. 


Similarly to July 2013 in June 2014, crabs were 
scanty and the sample was even smaller than in the 
previous year. The analysis has shown three size groups 
with mean CW 16.1, 19.9 and 23.0 mm respectively 
(Fig. 8C). These sizes are generally larger than respec- 
tive sizes of the three groups in 2011 but the smallest 
group may be compared to the intermediate group of 
2011, and the next to the largest group of 2011. The 
smallest group was strongly dominating; the other two 
consisted of few large crabs. 

Characteristically sampling in shallow coastal flats 
near Sennoi revealed few small crabs and can charac- 
terize only two or three adult size groups which have 


been persistently present in the population between 
2011 and 2014. 

Smallest size group in the Taman Bay 

Owing to occasional collection of small juveniles 
by bottom samplers in the benthic survey in July 2008 
(see above) we can characterize the smallest size 
groups presumably right after settlement. This group 
was well separated from few specimens that apparent- 
ly belonged to an older and larger group which was 
also present in the samples. Their mean CW was 3.6 
(Fig. 8A). Three juvenile crabs with CW 2.4-3.2 mm 
have been collected from Zostera marina at depth 3— 
4 m using SCUBA diving in July 2011. We designate 
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Fig. 6. Results of the mixing distribution analysis of size (CW, mm) composition of Rhithropanopeus harrisii collected in and out of 
burrows near Sennoi (Taman Bay) in July 2011. (see Table 8 for details). 

Puc. 6. PesyibraTeI pasgenenna CMellaHHoro pactipeyeneHua pasMepos (CW, MM) Ha HOpMaJIbHO paciipeyeueHHble Tpynnbi JIA 
Rhithropanopeus harrisii, coOpanubix B n BHe Hop y Toc. CenHolt (TaMaHckuit 3ammB) B mone 2011 r. (aeram cm. Tam. 8). 


therefore juvenile crabs with CW 2—5 mm in summer 
to group 0. 

Estuarine population in river Vulan in June 2013 

Male size population structure analysis has resulted 
in 3 groups, with mean CW sizes of 8.1, 14.7 and 19.4 
mm. (Table 10; Fig. 9A). The total female sample in 
the river has shown 3 groups (mean CW 8.8, 13.3 and 
15.7 mm) while females with and without eggs on 
pleopods taken separately comprised of only two groups 
(Table 10; Figs 9B—D). Groups of smaller size always 
outnumbered larger crabs. 

Neither these data nor the size composition of small 
crabs in the benthic survey of 2008 can be directly 
compared to the Sennoi transect data because of differ- 
ent methods of collection. However, it appears that a 
group with average CW around 8—9 mm detected in 
widely dispersed samples in June—July likely repre- 
sents the age group which is one year older than 0- 
group. At least part of females can mature and repro- 
duce at the age of around one year (Figs 7D, 9D). The 
age of larger size groups must be greater but can’t be 
estimated with certainty while the largest males (CW > 
18 mm) may have the age of more than 2 years. 


Dead crabs occurrence 


Dead crabs occurrence observed in July 2011 near 
Sennoi was of particular interest since it indicated sig- 
nificant mortality in the population in a restricted time. 
Out of 8 transects, 7 accounted for found dead crabs. 
Overall 42 dead crabs have been found out of burrows 
(mean CW size 16.33): 35 males (83%), 5 (12%) non- 


ovigerous females, and 2 (5%) ovigerous females. This 
proportion was not much different from the one found 
in live population. 

Crabs were largely intact and not yet in process of 
decomposition that indicated their relatively recent death 
(not longer than a day earlier). Mean CW of dead 
males was shifted towards the largest size group in the 
male part of population, while female sizes correspond- 
ed to the largest size group of the female part of popu- 
lation (Figs 10A, B). Overall the size distribution of 
dead crabs was strongly asymmetric with mostly large 
crabs dying. Thus predominately large groups were 
affected by mortality. We speculate that this could be 
due to natural ageing process probably associated with 
the impact of unknown pathogens. 

The populations densities in the following years 
have been much smaller and no dead crabs have been 
found. This once again suggests that the dead speci- 
mens found in 2011 could be a result of natural mortal- 
ity which was detectable only at high general abun- 
dance of the population. 


Discussion 


Distribution of Harris mud crab in the Black 
and Azov seas 


Vector of introduction of R. harrisii to the Black 
and Azov seas is not known with certainty. It was 
supposed that it came to the area with ship fouling 
[Makarov, 1939; Reznichenko, 1967]. Although, this 
species was recorded in the north-western Black sea 
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Table 6. Number of Harris mud crabs found in and out of burrows on different standard transects at different time of day in 
July 2011, in the Taman Bay (100% by day and transect In + Out). 

TaOmuua 6. KonmuectBo kpados Xappuca, OOHapy2KeHHEIX B HopKax (In) n BHe HOpoK (Out), B pa3Hoe BpeMa CyTOK Ha 
CTaHJapTHBIxX TpaHcekTax B urose 2011 r., B TamaHcKom sambe (100% ropH30HTaIbHO IMO HAM H TpaHceKkTaM B +BHe). 


Counts Morning Noon Evening 
Transect In Out In Out In Out 

1 101 (74%) 35 (26%) 94 (78%) 27 (22%) 73 (72%) 29 (28%) 
2 74 (40%) 113 (60%) 78 (55%) 65 (45%) 134 (77%) 41 (23%) 
3 4(41%) 66 (59%) 59 (49%) 62 (51%) 


more than a decade earlier than in the Northern Sea of 
Azov (1937 vs. 1949), we could not be sure that Harris 
mud crab had initially invaded the north-western Black 
Sea and only then entered the Sea of Azov. Both areas 
have large merchant harbors such as Odessa, Kherson 
and Nikolaev (the nearest port to the first recorded 
mud crab locality in the area) in the first area, and 
Mariupol, Taganrog and Rostov in the Taganrog Bay 
of the Sea of Azov. Introduction via shipping could 
also have happen at about the same time (i.e. in the 
1930s) in both areas, but the species was by chance 
first found at the advent of the World War II in the 
Dnieper-Bug liman and the war made research in the 
Sea of Azov impossible for several subsequent years. 

The northwestern Black Sea (from the Danube del- 
ta to the isthmus of the Crimean Peninsula), and the 
Sea of Azov are similar in regard to the Harris mud 
crab distribution pattern. There, this species occurs 
along the open continental shore, in offshore banks, in 
limans, lagoons and estuaries and ports [Mordukhai- 
Boltovskoy, 1952; Reznichenko, 1967; Makarov, 2004; 
Vinogradov et al., 2012]. In the Sea of Azov this spe- 
cies occurs in all habitats with the exception of oxygen 
depleted mud habitats in the deepest (10 m and deeper) 
area [Makarevich et al., 2000; Nabozhenko et al., 2010]. 
In some areas, i.e. the entrance to the Taganrog Bay, R. 
harrisii has been even recognized as a dominant spe- 
cies in the community according to the data of a small 
grab survey [Makarevich et al., 2000; Sergeeva, Bur- 
katsky, 2002]. 

In the Taman Bay our observations indicate similar 
pattern of Harris mud distribution throughout the bay, 
although the occurrence of crabs in its extremely shal- 
low north-western part [lacking sea grass vegetation; 
Spiridonov et al., 2016] is not yet confirmed. 

The first finding of the species in the northern part 
of the Kerch strait dates to 1955 [in an unpublished 
master thesis by Nesis, 1956]. By 1989 R. harrisii was 
a common species in the communities dominated by 
Mytilus galloprovincialis and Mya arenaria in the north- 
ern and central but not southern parts of the strait 
[Ivanov, Sinegub, 2008]. 

Beside our survey, several other detailed surveys of 
Crimean inlets failed to find this species on the Black 
Sea coast of this peninsula [Revkov et al., 2007; Alyo- 
mov, 2012] even though it was earlier reported from 
the inner part of Sevastopol Bay [Shalovenkov, 2005]. 


It appears that, along most of the Crimean and Cauca- 
sian coast of the Black Sea R. harrisii occurs patchi- 
ly. These patches are associated either with river estu- 
aries (Vulan, Shap’sukho) or lagoons (such as Bugas 
Liman), which have restricted water exchange with 
the sea. Larvae of R. harrisii were also reported from 
the harbors of Anapa, Novorossiysk and Tuapse [Seli- 
fonova, 2012]. The species has not been reported 
from Turkey’s coast of the Black Sea [Gönlügür- 
Demirci, 2006; Balkis et al., 2012]. No recent litera- 
ture sources on the decapods distribution are avail- 
able for Georgia’s coast. 

Even with significant knowledge gaps the distribu- 
tion of Harris mud crab in the Black and Azov seas 
appears to be highly uneven. The rarity of the species 
along the Crimean coasts resembles its known occur- 
rence pattern in the Mediterranean where R. harrisii is 
found in few lagoonal or estuarine habitats. 


Table 7. Two way ANOVA p results (factor “place” — in 
vs out of burrows, factor sample — samples, place : sam- 
ple — interaction term) of sizes of Harris mud crabs found 

on standard transects in the Taman Bay in July 2011. 
Tamna 7. PesynbTaTbI WByXPakTOpHOTO AHcrepcHOHHOrO 
anamu3a (1-1 daxtop “place” — B HJIH BHe HOpKH, 2-1 
tpaxtop “sample” — Homep mpo6sr, place:sample — 
B3aHMOelcTBHe (þakroposB) pa3sMepoB kpaðos Xappuca, 
coOpaHHBbIxX Ha CTaHJapTHBIxX TpaHcektax B TaMaHCKOM 
3ayIMBe B HHE 2011 r. 


Sex/group Factor Morning Noon Evening 
place 0.0797 0.92 0.784 
Males sample 0.0586 0.884 0.271 
place:sample 0.1142 0.63 0.49 
Non- place 0.00017 0.551 0.0661 
ovigerous sample 0.64443 0.537 0.4776 
femal 
emaes | place:sample 0.01678 0.859 0.2082 
place NA 0.8 0.9556 
Ovigerous. T sample NA 0.392 0.5585 
females 
place:sample NA 0.474 0.0953 


NA — Not applicable; * statistically significant; degrees of 
freedom is | in all cases except sample factor in the evening where 
df = 2. 

NA — HenlpHiioxKuMo; * cTaTHCTH4ecKH 3HAYHMOE BIIMAHHEe 
(paxtopa; KouMmuecTBO CTereHeli cBoOognI (df) Bo Bcex cly4asx paBHO 
1, kpome cbaxtopa sample B BeyepHee Bpema cyToK, rye df = 2. 
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Fig. 7. Results of the mixing distribution analysis of size (CW, mm) composition of Rhithropanopeus harrisii collected near Sennoi 


(Taman Bay) in July 2011 (see Table 9 for details). 


Puc. 7. PesyibraTsI pasgenenna CMellaHHoro paciipeyeneHua pasmepos (CW, MM) Ha HOpMaJIbHO paciipeyeueHHbIe Tpynnmbi JIA 
Rhithropanopeus harrisii, coOpanuerx y toc. Cenuoi (Tamaucxuii sanns) B mone 2011 r. (aeram cm. Taom. 9). 


The distribution pattern of the mud crab along the 
Black Sea coast of Romania and Bulgaria [Bacescu, 
1952; Petrescu, Balacsescu, 1995; Makarov, 2004; 
Skolka, Preda, 2010] and the Krasnodar area of Russia 
[present study] resembles the introduction area in the 
North Sea / Atlantic coast of Europe (Fig. 1A) and 
along the West coast of North America [Petersen, 2006], 
i.e. patchy association with estuaries and lagoons (or 
artificial lagoons such as harbors). 

The occurrence of R. harrisii throughout entire seas 
(i.e. Sea of Azov with the Kerch Strait), their large 


> 
w 


parts [north-eastern Black Sea, northern Caspian Sea 
— see Reznichenko, 1967] or big bays (Taman Bay) is 
not common for this species in other areas of introduc- 
tion. However this distribution pattern is found in its 
native range, along the eastern coast of North America; 
in particular in the Chesapeake Bay [Ryan, 1956] which 
is regarded as a large estuary. 

Significant part of the Sea of Azov is influenced by 
runoff of Don and Kuban rivers; the north-western 
Black Sea receives fresh water from three other big 
European rivers: Danube, Dnieper, and Dniester. Only 
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Fig. 8. Results of the mixing distribution analysis of size (CW, mm) composition of Rhithropanopeus harrisii collected during the 
benthic survey in the Taman Bay in 2008 (A) and near Sennoi town in 2013 (B) and 2014 (C). 

Puc. 8. PesymbraTsI pasgenenna CMellaHHoro pacnpenmenenna pasMepos (CW, MM) Ha HOpMaJIbHO paciipeyeueHHbIe Tpynmnmbi JIA 
Rhithropanopeus harrisii, coOpanubix B Xone OeHTOCHO! CÞeMKH B TaMaHcKom sanBe B 2008 r. (A) n y moc. Cennoii B 2013 (B) u 2014 


(C) rr. 
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Table 8. Characteristics of size groups revealed by the mixing distribution analysis of size composition of Rhithropan- 
opeus harrisii collected in and out of burrows near Sennoi (Taman Bay) in July 2011. F — Non ovigerous females, Fov — 
Ovigerous females, in — found in burrows, out — found out of burrows. 

TaOmuua 8. XapakTepHcTHKM pa3MepHBIX rpynl KpaOos, BbIABJICHHbIC C TOMOMI[bIO pa3esIeHiA CMelaHHOro 
pacipeyemenua pasmepos Rhithrpanopeus harrisii, coOpanuerx okozo Toc. CenHoli (TaMancknii 3aIMB) B HOpKax H BHe 
HOpOK B Hoe 201 Lr. F — camxu 6e3 nkphi, Fov — cCaMKH C UKpoii, in — B HOpKax, out — BHe HOpOK. 


Taman Bay 2011 Group Means (CW, mm) 
Log 

Sex/group N of groups likelihood Akaike IC 1 2 3 4 5 
Males in 5 —642.8 1306 7.6+0.6 9.6+0.5 11.7+0.6 15.2+1.6 18.6+1 
Males out 2 —308.6 625.4 12.1+1.3 15.741.6 
Fin 3 44.4 102.8 9.1+0.7 11.340.6 13.9+0.5 
F out 1 53.27 110.7 11.941.4 
Fov in 2 —50.57 109.8 9.5+0.4 12.241.1 
Fov out 1 -7.289 19.67 12.11 


Table 9. Characteristics of size groups of Harris mud crab revealed by the mixing distribution analysis of size composition 
of Rhithropanopeus harrisii collected near Sennoi (Taman Bay) in July 2011 (total sample). 

Tamua 9. XapakTepucTAKH pa3MepHBIX rpyn KpaOos, BÞIABJICHHDIE C MOMONDIO pa3esIeHiA CMEMAHHOTO 
pacnpegenennsa pasmepos Rhithropanopeus harrisii, coOpanutx okoo moc. Cennoŭ (Tamanckuň samb) B none 2011r. 
(o6pemunennas BEIOOpKa). 


Taman Bay 2011 Group Means (CW, mm) 

Sex/group N of groups Log likelihood_| Akaike IC 1 2 3 4 5 
All males 3 —962 1936 12.4+1.8 15.6£1.4 | 18.540.9 
All females 2 -173.1 354.4 9.2+0.7 12.11.2 
AIF 2 —107.9 224.3 8.80.5 11.9414 
All Fov 2 —58.8 126.1 9.6+0.4 12.2+1.1 


F — non-ovigerous females, Fov — ovigerous females. 
F — camxu 6e3 ukppi, Fov — camku C uKpoii. 


Table 10. Characteristics of size groups revealed by the mixing distribution analysis of size composition of Rhithropan- 
opeus harrisii collected from the river Vulan estuary (near village Arkhipo-Osipovka) in July 2013. 

Tamua 10. XapakTepucTHKH pa3MepHBIX Ipyliil, BbIABJICHHbIe C MOMOLIbIO pa3{eJIeHHA CMeLIaHHBIX paciipeyeseHnit 
pasmepos Rhithrpanopeus harrisii, coOpanubix B 9cTyapHu peke Byxan (B paiione nocemka Apxuno-OcnmoBka) B Hose 


2013 r. 
River Vulan 2013 Group Means (CW, mm) 
Sample size 

Sex/group N of groups Log Ik.hood | Akaike IC 1 2 3 

All Males 3 —61.9 138.9 8.141.9 14.742.1 19.4+0.4 34 
All Females 3 —46.7 109.3 8.842 13.340.2 15.7+41.9 29 
F 2 -20.2 52.9 6.9+0.9 13.342.2 14 
Fov 2 —18.6 49.25 9.740.4 1442.5 15 


F — Non-ovigerous females, Foy — ovigerous females. 
F — camxu 6e3 ukppi, Fov — camku C uKpoii. 
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the Taman Bay lacks river discharge. Although neither 
the Sea of Azov, nor the north-western Black Sea can 
in strict sense be regarded as very large estuaries, and 
the Taman Bay is not an estuary, because it lacks river 
runoff and salinity gradient, the characteristic species 
in benthic communities in these waters are typical for 
the transitional (estuarine-lagoonal) European habitats 
[Vorobyev, 1949; Nesis, 1956; Golovkina, Nabozhen- 
ko, 2012; Spiridonov et al., 2016]. Both in the shallow 
north-western Black Sea and the Sea of Azov currents 
are generally wind-driven and change direction in re- 
sponse to meteorological situation [Zalogin, Kosarev, 
1999]. Owing to this wind-driven circulation in semi- 
closed basins larvae of R. harrisii are maintained with- 
in a vast shallow area where salinity range does not 
exceed 18%o (and usually ranges between 10 and14 %o) 
and depth and substrate conditions may favor their 
settlement. In the Taman Bay, post-metamorphic young 
crabs can (as we found in 2008: Fig. 1 C) settle in 
communities dominated by bivalves and seagrass that 
provide shelter; it is also similar for the Sea of Azov 
and the northern Kerch Strait where in particular cock- 
le (Cerastoderma) and mussel (Mytilaster and Mytilus) 
beds are common [Vorobyev, 1949; Nesis, 1956; Iva- 
nov, Sinegub, 2008]. Furthermore the Sea of Azov and 
the Taman Bay are practically devoid of other crabs 
species which can be active predators/competitors for 
mud crabs (only Pilumnus hirtellus occurs in the Taman 
Bay). 

In coastal areas of the Black Sea along most part of 
the Crimean Peninsula and the Caucasian (north-east- 
ern and eastern) coast the shelf is narrow, while cur- 
rents usually have a dominant direction towards the 
north-west [Zalogin, Kosarev, 1999]. No larval reten- 
tion mechanism [Forward, 2009] which would act with 
sufficient predictability, can be supposed for these non- 
tidal areas, except for partial isolation of water bodies 
themselves (i.e. Bugas Liman) or existence of second- 
ary circulations in slow flow estuaries with islands (i.e. 
Shap’sukho and Vulan rivers). The larvae that occur in 
the waters near the open Crimean and Caucasian coasts 
are subjected to the risk of being transported to areas 
with greater depth (over 20-30 m) which can be subop- 
timal for settling of metamophizing larvae. 

In broad shallow area adjacent to the entrance of 
Kerch Strait neither depth nor salinity (up to 18—19%bo) 
are essentially different from the Taman Bay. However 
R. harrisii has not been found there (southern side of 
Tuzla spit, Zheleznyi Rog cape, off Blagoveschenskaya 
village and near Anapa; Fig. 1B) even though it occurs 
on the northern side of Tuzla spit and in the Bugas 
Liman which have direct connection to this part of the 
Black Sea (in the latter even at higher salinity than in 
the neighboring sea). Furthermore larvae of Harris mud 
crab have been found in the harbor of Anapa [Selifono- 
va, 2012]. We speculate that either greater exposure of 
the coast to wave action and the dominance of sandy 
sediments in the nearshore zone from Anapa to the tip 
of Tuzla spit makes this habitat unsuitable for R. har- 


risii or /and the presence of other larger crab species 
(especially Liocarcinus vernalis and Eriphia verru- 
cosa) increases predation risk for juvenile mud crabs. 
In fact, the Sea of Azov shores of Kerch Peninsula, 
where we have observed numerous mud crabs (Figs 
1B, 2A) are also open to high energy wave activity and 
the sediments are often disturbed, hence they are most- 
ly sand or shell debris, similar to the Black Sea shore 
along the Bugas Liman. However, in the Sea of Azov 
the juveniles may settle and live elsewhere, where the 
conditions are more suitable, while in the Black Sea 
adjacent to the Kerch Strait such conditions may be 
absent. It is very likely that the presence of sheltered 
habitats for settled juveniles, such as sea grasses in 
Taman Bay, Mytilaster beds in the Sea of Azov or 
oyster reef in native habitats [Tolley et al., 2013] are of 
primary importance for maintaining mud crab popula- 
tions in coastal areas. 

In conclusion to the discussion of R. harrisii distri- 
bution in the Black Sea let us note that while estuarine 
and lagoonal populations along the Romanian and Bul- 
garian coasts could originate either owing to larval 
transport with south-westward currents from the north- 
western Black Sea [Zalogin, Kosarev, 1999] or due to 
the shipping vector, i.e. vessels called to such ports as 
Constantza and Varna, the origin of Harris mud crab 
populations discovered by us in the estuaries of rivers 
Vulan and Shap’sukho is unclear. It is unlikely that the 
founders of these populations were brought by ships 
because there are no harbors in the vicinity of their 
localities. We suppose that the initial mud crab coloni- 
zation of estuaries and limans along the North-Cauca- 
sian coast of the Black Sea was due to larval transport 
by currents. Currents there are generally directed north- 
westward towards the Kerch Strait [Zalogin, Kosarev, 
1999] that makes advection of larvae from the Sea of 
Azov and the Taman Bay to the Caucasian coast of the 
Black Sea problematic. Alternative explanation would 
be an existence of yet not discovered populations of R. 
harrisii in the eastern Black Sea. It is likely that popu- 
lations of Harris mud crabs exist in protected harbors 
of Tuapse and Novorossiysk, because the larvae of this 
species have been found there [Selifonova, 2012]. These 
areas along with the harbor of Poti in Georgia could 
host a not yet discovered stock of mud crabs which 
might have provided propagule pressure sufficient to 
establish smaller populations in the estuaries of the 
north-eastern coast of the Black Sea. To test our hy- 
pothesis one needs to survey harbors of the eastern 
Black Sea coast and obtain molecular genetic data to 
analyze relationships between different Black/Azov seas 
populations of R. harrisii. 


Habitats and hiding habit 


Although known by common name as “mud crab” 
R. harrisii lives in variety of habitats. In the Azov and 
Black seas they include practically all existing types of 
substrate from muddy sediments to sand with shell 
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Fig. 9. Results of the mixing distribution analysis of size (CW, mm) composition of Rhithropanopeus harrisii collected from the river 


Vulan estuary in July 2013. (see Table 10) 


Puc. 9. PesymbraTeI pasgenenna CMellaHHoro pactipeyeneHua pasMepos (CW, MM) Ha HOpMaJIbHO paciipeyeueHHble Tpynnbi JIA 
Rhithropanopeus harrisii, coOpanuetx B 9cTyapuu peku Bynan B mone 2013 r. (cm. Taos. 10). 


debris and stones (in particular on the Sea of Azov 
coast of Crimean Peninsula near Zolotoe village, see 
Fig. 2A). R. harrisii is spread throughout the Taman 
Bay and is found in all biotopes. At maximum depths 
(over 3 to 5 m) crabs live on muddy sediments with 
shell. Females and juveniles were observed hiding un- 
der shells or in the folding of bacterial biofilms. In 
Zostera and Potamogeton overgrowths the crabs are 
hiding in burrows that are most probably dug out by 
gobies. Females are found in small inconspicuous pock- 
ets created by macrophyte roots with no definite en- 
trances. At no time have we observed the crabs creat- 
ing these burrows. The only burrowing activity that has 
been observed is a quick burrowing in the sand when 
attempting to hide from the observer. It seems that the 
crabs use available crevices and burrows created by 
other organisms or natural formations as in clay mazes 
or under stones and shells. In rivers, the crabs use roots 
of plants and stones that are available along the shore. 
It is unclear what happens on the bottom due to very 
low visibility, but the bottom of Vulan river is very soft 
and releases gas when disturbed. It is unlikely that 
there are any burrows in such sediments. 

Our observations are in agreement with reports of 
other authors that studied Harris mud crab in both 
native and introduction areas. It has been reported that 
crabs live in oysters reefs (in native habitats), under 
stones, in various plant debris, in mud and sand with 


shell, and even among Fucus macroalgae on hard bot- 
tom [Ryan, 1956; Turoboyski, 1978; Petersen, 2006; 
Fowler et al., 2013]. In the experiment they preferred 
vegetated boulders compared to unvegetated boulders 
or sandy habitats [Nurkse et al., 2015]. This indicates a 
highly opportunistic habitat choice; however the pres- 
ence of shelter appears to be a key factor [Petersen, 
2006]. 

Our observations indicate that the use of shelter and 
hiding pattern significantly vary between different size, 
sex groups and localities. In the shallow sandy and 
shelly coastal beds with macrophyte overgrowth of the 
Taman Bay considerable part of population remains in 
the open while males and females showed different 
pattern of hiding (Tables 5, 6, 7 and Figs 4—6). Majori- 
ty of males hide during the day and there is a tendency 
for smaller males to sit in burrows more often. Unfor- 
tunately, we were unable to conduct comparable night 
counts of crabs in the shelter and in the open. However, 
almost equal proportion of males in and out of burrows 
in the morning indicate greater occurrence of the crabs 
in the open in the dusk and/ or dark time which is 
probably associated with foraging. Contrary to this, 
non-ovigerous females preferred to be out during the 
day and hide towards the evening. In the morning, 
when significant part of males was still outside, larger 
females were out of shelters and smaller were hiding in 
the burrows. This indicates a separation of peaks of 
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Table 11. Maximum and minimum sizes (CW, mm) of all Harris mud crabs crabs, males and females found on standard 
transects in the Taman bay at different years and in the river Vulan in 2013. 

TaOmua 11. MakcuMasIbHbie H MMHMMAaJIbHBIe pa3sMepbl (LIMpHHa Kapamiakca, MM) Bcex KpaðoB Xappuca, CaMI{OB H CaMOK, 
OOHapy2KeHHBIX Ha CTaHapTHbIX TpaHceKTax TaMaHCKOrOo 3aJIMBa B pa3Hble TOMBI, a TAKKE B ICTYapHU pekun Bynas B 2013 


P 
Locality | Yem | crate | M®xCW | Cw | males | males. | females | females 
Jul 2011 1096 20 7 20 7 15.5 8 
Taman Bay Jul 2013 64 23.9 10 
Jun 2014 21 23 13.1 
River Vulan July 2013 63 20 5 20 5 19.1 5 


Samples in the Taman Bay in 2013 and 2014 consisted mostly of males. 
Ilpoő6s: n3 Tamaucxoro 3asuBa B 2013 u 2014 rr. COCTOAJIM npeHMyllecTBeHHO H3 CAMIOB. 


outside activity in crabs of different sex that probably 
leads to a decreased frequency of agonistic encounters 
between females and males which generally have larg- 
er size. Ovigerous females generally had cryptic habit 
during day light. Even though their out of burrows 
sample was too small for statistical analysis, one could 
speculate that while smaller ovigerous females were 
still hiding in the morning, the larger ones showed 
greater activity in open space. The mesocosm experi- 
ments with Harris mud crabs in the Baltic showed that 
they also have size-associated pattern of using differ- 
ent types of substrates and shelter and interspecific 
aggression and food availability play a role in crabs’ 
pattern of using different habitats [Nurkse et al., 2015]. 

In contrast to R. harrisii in the Taman bay crabs in 
the Sea of Azov (near Zolotoe village) and in the river 
estuaries were practically not observed in the open. 
Crabs in the estuaries hide whenever there is slight 
movement. Such behavioral difference may be due to 
differences in predators’ pressure in studied areas. Not 
much is known about predators of mud crabs in the 
Azov and Black seas except that they are consumed by 
several fish species [Reznichenko, 1967]. According 
to our visual observations potential predators in the 
shallow areas of the Taman Bay are mainly gobies usu- 
ally not exceeding 30 cm in length and occur in moder- 
ate abundance (not more than 2—5 specimens per 10 m°). 
Cases of sharing burrows between gobies and crabs may 
indicate that at least for large adult crabs gobies are not a 
serious danger and combination of hiding (using exist- 
ing burrows, macrophyte shelter and rapid burrowing in 
sand) and taking a defense position may be an effective 
strategy to tackle with potential predators. 

Near Zolotoe village only stones provide shelter 
while sediments (shell debris and coarse sand) may be 
hard for rapid burrowing. Together with high abun- 
dance of large size gobies (Fig. 2A) this probably 
resulted in a totally adopted cryptic habit of crabs. 

We can only speculate that the hiding behavior of 
Harris mud crabs in the north-eastern Black Sea estuar- 
ies is also related to the interaction with predators. 
According to communication of local people freshwa- 
ter fishes are common in the estuaries thus possibly 
provide high predator pressure. Other potential preda- 


tors could be water birds. Not only the ability to use a 
variety of substrates, but also a flexible hiding behav- 
ior is part of opportunistic strategy of R. harrisii facili- 
tating its survival and thriving in various areas of intro- 
duction. 


Size range of Harris mud crabs 


Size composition data of R. harrisii in literature 
came from using different methods of collecting, so it 
is only the size ranges that can be generally compared. 
In different localities of Poland’s Baltic waters the 
minimum size (CW) of males ranged between 2.9 and 
5.6 mm and the maximum size — between 21.4 and 
26.0; those of females ranged between 1.9-5.3 and 
16.1-19.8 mm [Hegele-Drywa et al., 2014]. In the Bay 
of Finland, respective ranges for males and females 
were 2.2—22.0 and 2.1—20.2 [Fowler et al., 2013].The 
size range of crabs in the Taman Bay and in the Vulan 
River estuary (Table 11) corresponds well with the 
data for the Baltic Sea populations. It was earlier ob- 
served [Fowler et al., 2013; Hegele-Drywa et al., 2014] 
that R. harrisii attains greater size in its Baltic intro- 
duction area than in its native distribution range, i.e. 
Chesapeake Bay, Louisiana where males larger than 18 
mm, and females larger than 16 mm have not been 
reported [Ryan, 1956; Williams, 1984; Fowler et al., 
2013]. In our Taman Bay samples CW around 18 mm 
was repeatedly recorded as a mean size of the largest 
size group (Figs 6—8; Tables 8, 9). On the basis of our 
data and observations by Makarov [1939] we conclude 
that attaining larger size in introduction areas also stands 
true for the Black and Azov seas population of Harris 
crab. It follows a general pattern revealed for most 
studied marine invasive species, the majority of which 
become significantly larger in the introduced range 
compared to the native range with little evidence for any 
decrease in size following invasion [Grosholz, Ruiz, 
2003]. Even though no mechanism explaining this in- 
crease has been demonstrated in details [Grosholz, Ruiz, 
2003] our observations in the Taman Bay indirectly 
suggest that the decrease of predator pressure may be an 
important factor. This should be taken into account in 
further studies of mud crabs in the Black and Azov seas. 
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Fig. 10. Box plot of size versus sex (A) and size histogram (B) of dead crabs found on standard transects in the Taman Bay in July 
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Puc. 10. J[uarpamma pa3maxa pa3Mepos (A) m rucrorpamma (B) pa3MepHoro pasga MepTBBIX KpaOoB (pa3evIeHHBIX IMO MOJY), 
coOpaHHBIX Ha CTaHapTHBIX TpaHcekTax B Hose 2011 r. B TamaHcKoM 3aJIuBe. 


Composition and dynamics of Harris mud crab 
in coastal shallow beds of the Taman Bay 


The shallow coastal flats near Sennoi where we 
conducted transect censuses of R. harrisii in 2011, 
2013 and 2014 appear to be generally favorable habi- 
tats for this species in summer. This is due to relatively 
low predators pressure, at least for larger crabs (see 
above) and a mosaic pattern of macrophyte dominated 
benthic communities providing a variety of prey items 
for Harris mud crabs which are generalist predators and 
scavengers [Zalota et al., 2013; Kolesnichenko, 2014]. 
The disadvantages for living in this zone include risks of 
storm disturbance, particularly in autumn and winter, 
and heating of the water in calm summer days which 
may cause hypoxia recorded from time to time in shal- 
low area of the Taman Bay [Spiridonov et al., 2016]. 

Male to female ratio in our collections in the Taman 
Bay shallows was biased towards males 78%; Table 4) 
in 2011 while in 2013 and 2014 samples consisted 
practically entirely of males. Broad surveys in Dead 
Vistula and the Gulf of Gdansk using dredges indicated 
only slight excess of males. i.e 54% [Turoboyski, 1973] 
and 52% [Hegele-Drywa et al., 2014] of adult popula- 
tion respectively. Similar fraction of males has been 
reported by Reznichenko [1967] for the Sea of Azov in 
general (although the method of collecting was not 
described in details). In the season-long crate collect- 
ing in the Gulf of Finland male to female ratio varied 
from 44 to 51% and did not show significant differenc- 
es between locations and different years [Fowler et al., 
2013]. On the other hand in the Taganrog Bay males 
comprised 82% (similar to our results) of samples 
[Reznichenko, 1967]. Another surprisingly similar fig- 
ure (about 82% of males) was found by Kolesnichenko 


[2014] who collected crabs by digging sediments near 
the watermark in the Russian part of the Vistula Bay in 
summer. Vorstman [1939] also reported several cases 
of strong (up 100%) dominance of males in particular 
habitats in Netherlands. 

One of the explanations for biased sex ratio in the 
coastal shallows of the Taman Bay is that females 
which maintain cryptic way of life were underestimat- 
ed. In this case we would expect a much higher number 
of females found in the shelter (burrows) than out of 
the shelter. This was not however the case (Table 5), so 
the main source of underestimation could be that the 
females completely burrowed themselves in sediments. 
However if this was a common case we would observe 
at least certain amount of crabs sitting entirely in sand 
because the inspection of transects proceeded with slow 
speed and was carefully surveyed. Indeed some crabs 
burrowed themselves in sand when being disturbed but 
this was easily noticed. Therefore we believe that under- 
estimation of females was not really high and our figures 
reflect a real situation in the shallows of the Taman Bay. 

Our study revealed two or three adult size groups of 
Harris mud crab which have been persistently present 
in the coastal shallows of the Taman Bay between 
2011 and 2014. Characteristically smaller size groups 
were deficient. This could again be due to the collec- 
tion method but in 2011 small crabs of 5-6 mm (CW) 
were spotted and collected. We could also see and 
potentially count other small creatures such as hermit 
crabs, /dothea, and others. Thus even assuming under- 
estimation of smaller size groups of crabs we conclude 
that shallow coastal areas in the Taman Bay in summer 
are populated by predominately males of larger size 
groups. Females may be also present as we discovered 
them in summer 2011, some of them being in process 
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of carrying developing eggs on pleopods, and in April 
2014 when mating was observed [unpublished obser- 
vations of the first author]. However, one can easily 
suggest that release of larvae in the shallow banks 
poses high risks of negative effects of surge on larvae. 
Thus larvae hatching in this part of the bay could be an 
unsuccessful strategy and females most likely migrate 
to deeper areas prior to hatching. Consequently during 
the breeding season sex ratio of crabs inhabiting coast- 
al shallows is biased towards males. Also coastal flats 
near Sennoi (depth < 1.5 m) are unlikely to be populat- 
ed by crabs in winter because of frequent storms and 
seasonal diminishing of macrophyte vegetation [obser- 
vations by authors]. Therefore, it is likely that the pop- 
ulation of R. harrisii in the Taman Bay has a complex 
spatial structure with different age groups, males and 
females occupying partly different habitats and per- 
forming seasonal migrations. Even though our sam- 
pling methods in the Taman Bay and in the river estu- 
aries do not allow for full correct comparison of the 
data, simultaneous presence of juvenile crabs and adult 
females (including ovigerous) and males (Fig. 9) indi- 
cate that in the estuaries all these groups live more or 
less together. This is because the area populated by 
mud crabs is simply much too small for the establish- 
ment of a complex spatial structure. We can further 
speculate that spatial overlapping of different age groups 
and sexes in the estuaries result in frequent agonistic 
interactions between crabs and even cannibalistic pre- 
dation of smaller and molting specimens so that ex- 
tremely expressed hiding pattern of the estuarine crabs 
receive additional explanation. 

Our four years transect censuses in the Taman Bay 
covered only part of population which temporarily oc- 
cupies coastal shallow bank. The timing of occupation 
is also not known with certainty. Only with caution we 
can interpret tremendous changes in abundance be- 
tween 2011 on one hand and 2013 and 2014 (by two 
orders of magnitude) on the other hand as an indication 
of pulsation of population abundance between years. 
When having frequently observed dead crabs in 201 lwe 
probably witnessed ageing and dying off a highly abun- 
dant generation, which exceeded in strength genera- 
tions of subsequent years. 


Further development of methodology to study 
Harris mud crab populations 


The method of direct observation and hand collec- 
tion of Harris mud crab in situ used in present study 
clearly shows its strength and limitations. First of all it 
is probably optimal for rapid assessment and detection 
of the species in relatively clear water with various, 
including coarse substrates and abundant macrophyte 
growth. On the coasts of the Black and Azov seas it is 
logistically most feasible because everything that is 
required is a vehicle and snorkeling/ diving equipment. 
It also provides valuable information on habitats and 
especially microhabitats populated by crabs and pat- 
tern of using shelter that would be practically impossi- 


ble to obtain by any other means, including dredging 
[Turoboyski, 1973] and long term collecting in crates 
[Roche, Torchin, 2007]. Underwater transect censuses 
also provides quantitative estimates of abundance but 
only for particular habitats, i.e. shallow areas with mixed 
sediment bottom and appropriate visibility. However, 
besides visibility there are several problems in applying 
transect counting, in particular massive development of 
filamentous algae on the bottom that prevents effective 
search of crabs. In deeper areas if SCUBA is used addi- 
tional technical problems arise from unavoidable dis- 
turbing of sediments that also worsen search condition. 
Using a less invasive and rapidly applicable to larger 
surveyed areas methods such as applying a calibrated 
bar for virtual construction of a sampling quadrate [Mil- 
jutin, 2005] can probably solve this problem. 

To study composition and dynamics of such exten- 
sively distributed and spatially structured populations 
of R. harrisii as the population in the Taman Bay a 
combination of methods should be applied with direct 
censuses in shallow areas, dredging on soft bottom 
with sparse vegetation, and qualitative collecting in 
dense Zostera marina beds using SCUBA. For estuar- 
ies similar to Vulan river, as well as for harbors the 
most appropriate way to collect material comparable 
between sites and years is installation of crates. Popu- 
lations of R. harrisii in other areas may require specific 
combinations of different sampling approaches. 


Conclusions 


The present study provided new data on the distri- 
bution of Rhithropanopeus harrisii in the eastern part 
of the Black Sea — Azov Basin and indicated its patchy 
distribution along the eastern coast of the Black Sea 
with newly discovered populations in two river estuar- 
ies (Vulan and Shap’sukho), and in the Bugas liman 
adjacent to the Black Sea. Furthermore we have shown 
a wide occurrence of the species in the Taman Bay 
located at the boundary between the Sea of Azov and 
the Black Sea and opening to the Kerch Strait. In the 
shallow area of the Black Sea adjacent to the Taman Bay 
and the Bugas Liman (where neither depth nor salinity 
are essentially different from the above areas) no indica- 
tions of Harris mud crab’s presence were found. We 
speculated that this may be due to limitation of optimal 
possibilities for settlement of young crabs and/ or the 
presence of abundant native crab species populations. 

Harris mud crabs have been found in very distinct 
shallow water habitats spanning from exposed shores 
with sandy, shell and stones on the Crimean coast of 
the Sea of Azov to extensive Zostera spp. meadows in 
the Taman Bay and narrow bands of reed roots along 
the river estuaries. The only common characteristic of 
these habitats is their ability to provide shelter. Direct 
observations and censuses in the Taman Bay indicated 
that smaller crabs of both sexes and most of ovigerous 
females prefer to hide in burrows. Overall, a flexible 
hiding behavior is part of opportunistic strategy of R. 
harrisii facilitating its survival and thriving in various 
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areas of introduction. Seemingly, there is a seasonal 
migration pattern of R. harrisii in the shallows of the 
Taman Bay and no such observations are apparent in 
the estuary habitats on the Black Sea coast. 

Similarly to the Baltic, R. harrisii attains greater 
size in the Black Sea and the Sea of Azov than in the 
native distribution range. A complex size structure in- 
cluding two or three size groups of adult/ subadult 
crabs (except newly settled juveniles) is characteristic 
for both the Taman Bay and the estuarine populations. 
Thus, although, the distribution and migration patterns 
of different size groups differ, the general life cycle of 
the studied crab populations appears to be similar. 
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